Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 2 February 2008 



(MN WT&L style file v2.2) 



XMM-Newton RGS spectroscopy of LMC X-3 



M.J. Page 1 , R. Soria 1 , K. Wu\ K.O. Mason 1 , F.A. Cordova 2 and W.C. Priedhorsky 3 

1 Milliard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey, RH5 6NT, UK 
2 Department of Physics, University of California, Riverside, CA93106, USA 
3 Los Alamos National Laboratory, Los Alamos, NM 87545, USA 



2 February 2008 



ABSTRACT 

We present soft X-ray spectroscopy of the black hole binary LMC X-3 from the 
XMM-Newton Reflection Grating Spectrometer. The observations span the full 
range of spectral states seen in LMC X-3. The spectra are completely dominated 
by continuum emission, and the neutral absorbing column measured from the 
O I edge (Nh = 3.8lg"7 x 10 20 cm" 2 ) is consistent with the Galactic interstellar 
column density towards LMC X-3. We find no evidence for variability of the neutral 
absorbing column. We also constrain the ionized column density using the upper 
limits to the equivalent width of the O II - O VIII K-shcll resonance lines: we 
find that the equivalent hydrogen column density of gas in which O is partially 
ionized is < 8 x 10 20 cm" 2 . From this upper limit we can rule out a line driven 
stellar wind as the power source for the X-ray emission of LMC X-3 except when 
it is faint. At wavelengths longward of the peak emission the spectral shape is 
well described by a multi-temperature disk-blackbody spectrum; the powerlaw 
component which dominates at shorter wavelengths does not continue longward of 
the disk-blackbody peak. This implies that the multi-temperature disk-blackbody 
component supplies the seed photons which are Compton upscattered in the hot 
corona, consistent with the standard paradigm for black hole X-ray binary spectral 
states. 

Key words: accretion, accretion disks - black hole physics - galaxies: Seyfert - 
binaries: general - stars: early type - stars: winds 



1 INTRODUCTION 

The X-ray properties of accreting black hole candidates 
(BHCs) lead naturally to classification into different spec- 
tral states (e.g. Van der Klis 1994). In their "high/soft" 
state, their X-ray spectra are dominated below 10 keV by a 
thermal blackbody-like component with temperature ~ 1 
keV. In their "low/hard" state, this thermal component is 
absent, and BHCs have powerlaw X-ray spectra. An "in- 
termediate" state is also seen, in which the thermal com- 
ponent is weaker. Individual BHCs can exhibit all spectral 
states, but the fraction of time spent in the different states 
varies considerably between objects: for example LMC X- 
1 has only been observed in the high/soft state, while Cyg 
X-l spends the majority of its time in the low/hard state. 

Of all the persistent BHCs known, LMC X-3 is viewed 
through the smallest Galactic column density. Until re- 
cently LMC X-3 was thought to be permanently in a 
high/soft state, but observations with the Rossi X-ray 
Timing Explorer (RXTE) have shown that LMC X-3 
changes to a low/hard state during periods of low count 
rate (Wilms et al. 2001). The existence of state transitions 



in LMC X-3 has now been confirmed by further RXTE 
observations (Boyd et al. 2000) and XMM-Newton obser- 
vations (Wu et al. 2001). Wilms et al. (2001) argue that 
the long term variability and consequent state changes 
could be due to an accretion disk wind-driven limit cycle 
(Shields et al. 1986). 

Because of the small column density in the direction 
of LMC X-3, it is the best persistent BHC in which to 
study the soft X-ray spectrum at a variety of spectral 
states. The soft X-ray spectral region is particularly im- 
portant for understanding the accretion and mass transfer 
processes, because winds from the accretion disk, or from 
the companion star will imprint the soft X-ray spectrum 
with emission and absorption lines and edges. Further- 
more, the standard paradigm for the thermal component 
seen in the high and intermediate states of BHCs (namely 
that it is the multi-temperature blackbody spectrum of 
an optically thick, gemoetrically thin, accretion disk) has 
been poorly tested at soft X-ray energies. This is because 
almost all instruments on the previous generation of X- 
ray observatories (i.e. before the launches of Chandra and 
XMM-Newton) had poor spectral resolution below 2 keV, 
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Figure 1. The 3-day averaged RXTE All-Sky Monitor 

lightcurve for LMC X-3 for the year 2000, with the times of 
the XMM-Newton observations marked. 



Table 1. XMM observations of LMC X-3 



Date 


XMM 


Instrument 


Exposure* 


Count rate* 


(2000) 


Orbit 




(ks) 


(count s — 1 ) 


07 Feb 


0030 


RGS1 


12.1 


7.1 


08 Mar 


0045 


RGS1+2 


20.8 


5.9 


19 Apr 


0066 


RGS1+2 


44.3 


0.013 


09 Jun 


0092 


RGS1+2 


77.5 


2.2 


24 Nov 


0176 


RGS1+2 


20.9 


18.0 



* Where applicable this is the mean of RGS1 and RGS2 

leading to degeneracies between different spectral shapes 
and different amounts of absorption in spectral modelling. 

In this paper we present soft X-ray spectra of LMC X- 
3 from the XMM-Newton Reflection Grating Spectrometer 
(RGS). The spectra, which have a resolution of ~70 mA, 
cover the full range of observed spectral states in this 
source. 



2 OBSERVATIONS AND DATA REDUCTION 

Table □ lists the XMM-Newton observations of LMC X-3 
that will be used in this analysis, and in Fig. Q the ob- 
servation times are marked on the RXTE All-Sky Mon- 
itor lightcurve for LMC X-3 during 2000. LMC X-3 was 
also observed by XMM-Newton in revolution 0028 on the 
2nd February 2000, but we have excluded this observa- 
tion from the present analysis because we are unable to 
produce a satisfactory effective area calibration for rev. 
0028, in which the source is off-axis. Over the ten months 
spanned by these observations LMC X-3 varied by more 
than a factor of 1000 in RGS count rate. 

The data were initially processed using the public re- 
lease XMM-Newton SAS version 5.2. First and second or- 
der spectra from both RGS were extracted (except for 
revolution 0030 in which only RGS1 was operating) and 
response matrices were generated using the latest version 
of the SAS task RGSRMFGEN (corresponding to SAS ver- 
sion 5.3.3). To correct for residual artifacts in the effective 
area calibration (see den Herder 2002), the effective area 
of each response matrix was then divided by the ratio of 
a powerlaw + Galactic column fit to the XMM-Newton 
rev 0084 RGS spectra of the continuum source Mrk 421 



(Brinkmann et al. 2001). LMC X-3 was on-axis for all of 
the observations listed in Tabled as was Mrk 421 in rev- 
olution 0084. The data/model ratio of the Mrk 421 RGS 
spectrum is shown in Fig. [5] for most of the RGS wave- 
length range the effective area correction is only a few 
percent. The statistical errors on effective area correction 
were added in quadrature to the statistical errors on the 
data. Finally, for each revolution, spectra and response 
matrices were resampled and coadded to produce a sin- 
gle spectrum per observation. To improve signal to noise, 
spectra from the February, March, June and November 
observations were grouped by a factor of 3, resulting in 
spectra with ~ 1000 channels ~ 30 mA wide, well sampled 
with respect to the RGS resolution of ~ 70 mA FWHM. 
The flux level of the April spectrum is extremely low, nec- 
cessitating very heavy binning before a respectable signal 
to noise ratio is achieved; we have grouped this spectrum 
by a factor of 64. The resultant spectra were then analysed 
using the XSPEC software package version 11.2.0 (Arnaud 
1996). 



3 RESULTS 

Fig. [3] shows the flux calibrated spectra from all five ob- 
servations. The four brightest spectra are shown at full 
resolution in Figs [I] -0 The spectra appear to be contin- 
uum dominated without prominent absorption or emission 
lines; this is particularly clear in the 24th November spec- 
trum, which has the highest signal to noise ratio and shows 
no significant spectral features except for the interstellar 
O I edge at 23A. 

3.1 X-ray absorption in LMC X-3 

The degeneracy between absorption and the intrinsic spec- 
tral shape has been a major source of uncertainty in at- 
tempts to model relatively low resolution soft X-ray spec- 
tra of BHCs (e.g. Haardt et al. 2001). The excellent combi- 
nation of sensitivity and resolution afforded by the XMM- 
Newton RGS allows individual absorption edges and lines 
to be measured more or less independently of the under- 
lying continuum shape, allowing a far more reliable deter- 
mination of the amount of soft X-ray absorption. 

3.1.1 Cold material 

Cold material (either interstellar or local to the source) 
imprints a well known sequence of absorption edges on 
X-ray spectra. For cosmic abundances, oxygen produces 
by far the strongest absorption edge in the RGS energy 
range, at 23 A, and inspection of Fig. [3] shows that this 
is consistently the strongest feature observed in the RGS 
spectra of LMC X-3. We have therefore used the depth 
of the 23 A O I K edge to determine the cold absorb- 
ing column for each of the spectra. By studying a small 
spectral region (3A to either side of the edge) we ensure 
that our column density determination is essentially inde- 
pendent of the overall continuum shape. The confidence 
intervals for the position of the edge and its optical depth 
are shown in Fig. [S]for all the observations excluding the 
very faint April observation. We have also fit the four 
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Figure 2. Data/model ratio for the RGS spectrum of Mrk 421 used to correct residual artifacts in the effective area calibration. 



observations simultaneously with a single set of edge pa- 
rameters but with the continuum slope and normalisation 
allowed to vary independently for each observation. The 
best fit value for the simultaneous fit is shown as a cross 
on the confidence intervals in Fig. |H| The edge parame- 
ters derived from each of the observations independently 
are consistent with the best fit edge parameters from the 
simultaneous fit; there is no significant change in the col- 
umn density of cold material from one observation to the 
next. The wavelength of the edge from the simultaneous fit 
is A = 22.98+007A where the confidence interval is 95%. 
This is consistent with the wavelength of the edge calcu- 
lated by Verner et al. (1996) , which is the edge wavelength 
used in the XSPEC 'tbabs' model (Wilms, Allen & McCray 
2000). To determine the equivalent H I column density we 
substituted the 'tbabs' model for the edge and performed 
another simultaneous fit to the four spectra. This results 
in Nh = 3.8^0 7 x 10 20 cm -2 , which is consistent with the 
H I column density in the direction of LMC X-3 as deter- 
mined from radio observations (3.2 x 10 20 cm -2 ; Nowak 
et al. 2001). 



3.1.2 Ionized material 

Although the column density of cold absorbing material 
is consistent with the interstellar medium of the Galaxy, 
absorbing material associated with a wind from the ac- 
cretion disk or the secondary star may be highly ionized 
by the intense X-ray and ultraviolet radiation. Absorption 
and emission lines from warm, ionized material have now 
been detected in the soft X-ray band from a number of 
AGN and X-ray binaries (e.g. Kaastra et al. 2000, Brandt 
& Schultz 2000). 

There are no obvious emission or absorption lines 
present in the RGS spectra of LMC X-3 (FigsgHHUniandQ 
and so we have attempted a sensitive search for their pres- 
ence. The interstellar medium and the Local-Group inter- 
galactic medium have imprinted absorption lines of O I at 
23.52A and O VII at 21.60A on the spectrum of Mrk 421 
(Paerels et al. 2003), which we have used to correct resid- 



ual artifacts in the effective area calibration. Therefore 
we performed the absorption line search on spectra which 
have no effective area correction. First, we calculated the 
transmission through the interstellar medium as a func- 
tion of energy using 'tbabs' with Nh = 3.8 x 10 20 cm -2 . 
We then corrected the spectra for interstellar absorption 
by dividing the effective area of the response matrix by 
the transmission. We then produced continuum spectra 
by smoothing our RGS spectra with a sliding box 30 pix- 
els (~ lA) wide. We rejected any points more than 2a 
deviant from the resulting smooth curve (i.e. potential 
emission/absorption features), and repeated the smooth- 
ing. This smoothing/rejection process was iterated several 
times to ensure that even marginally significant discrete 
features do not influence its shape. We then performed a 
grid search for emission and absorption lines where the 
grid was defined in terms of central wavelength, FWHM 
and equivalent width. The grid values of central wave- 
length were chosen to oversample the wavelength bins of 
the data, FWHM was stepped between and 1000 km 
s _1 , with a stepsize of 200 km s _1 (well sampled at the 
RGS resolution) , and the stepsize of the equivalent width 
was matched to the statistical uncertainty on the data. At 
each point in the grid, we formed a model spectrum by 
adding a gaussian emission/absorption line to the smooth 
continuum spectrum, and we computed the goodness of 
fit of this model to the real (unsmoothed) spectrum using 
the poisson errors on the data points and the x 2 statis- 
tic. The significance of any features, and limits on the 
equivalent width were determined using A^ 2 - Since we 
search for lines in either emission or absorption, we expect 
A% 2 to give an approximate measure of the significance of 
any lines (see Protassov et al. 2002). The grid search was 
performed using our own software which propogates the 
model through the response matrix in an identical fashion 
to XSPEC. 

The results of this line search are shown in Fig. |^] 
There are -~1000 channels in each spectrum, and there- 
fore we expect ~ 1 spurious feature at the 3er level in 
each spectrum. The spectra from the 7th February and 
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Figure 4. RGS spectrum from 07 February 2000 




Figure 5. RGS spectrum from 08 March 2000 



8th March each contain one 3a significant absorption fea- 
ture, at 36.18±0.08A and 32.61±0.05A respectively. These 
do not correspond to the wavelengths of strong tran- 
sitions in abundant elements, and are therefore proba- 
bly spurious. In the spectrum from 9th June, there are 
three 3a significant absorption features, at 17.86±0.08A, 
21.62±0.04A and 37.38±0.05A, and one 5a significant fea- 
ture at 23.51±0.03A. The features at 17.86A and 37.38Ado 



not correspond to the wavelengths of any expected ab- 
sorption line and are again probably spurious. The 21.62A 
feature has an equivalent width of 19±12mA and is consis- 
tent with the Is — 2p transition in O VII, while the 23.5lA 
feature is the Is — 2p transition of O I. In the spectrum 
from 24th November, the most significant feature (6a) is 
again the Is — 2p transition in O I. In addition, there 
are three 3a and one 4a significant absorption features 
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09 June 




Figure 6. RGS spectrum from 09 June 2000 




Figure 7. RGS spectrum from 24 November 2000 



in this spectrum. Two of the 3a features, at 14.94±0.04A 
and 31.30±0.06A do not correspond to strong transitions 
in abundant elements and are therefore probably spurious. 
The other 3cr line, at 18.93 ±0.05A, corresponds to O VIII 
Lya, while the 4cr line is found at 21.58 ± 0.03A, and cor- 
responds to O VII He-a. The 18.93A and 21.58A have 
equivalent widths of 9±4mA and 21 ± 12mA respectively. 



The O VII lines detected on the 9th June and 
24th November have similar equivalent widths of ~20mA, 
which implies an O VII column density between 
3xl0 15 cm -2 and 5xl0 16 cm~ 2 . Intervening material as- 
sociated with the halo of our Galaxy is likely to make a 
significant contribution to (and may dominate) this line, 
because Galactic O VII absorption is seen with a simi- 
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Figure 8. Confidence intervals for the depth and wavelength of the OIK edge for the February, March, June and November 
observations. The solid, dashed, and dotted lines correspond to 1,2 and 3<j respectively. The best fit edge parameters from a 
simultaneous fit to all four spectra is marked with a cross in each panel. 



lar equivalent width in the blazars observed during RGS 
calibration (Paerels et al. 2003). It is notable that if we 
apply the Mrk 421 effective area correction before per- 
forming the line search on LMC X-3, the OVII absorption 
lines are not detected in any of the spectra. The equivalent 
width of the O VIII Lya absorption observed on the 24th 
November implies an O VIII column density of between 
4xl0 15 cm -2 and 1.4xl0 16 cm -2 . This absorption line is 
seen in two out of three RGS calibration blazars (Paerels 
et al. 2003) with a similar equivalent width, and so the 
line in LMC X-3 may or may not arise in the Galactic 
halo. The observed wavelength is consistent with an out- 
flow from LMC X-3, but is not well enough constrained to 
rule out a Galactic origin. 

We can place very firm upper limits on the column 
density of an ionized wind by using the 3<r upper limits on 
the equivalent widths of Is — 2p absorption lines from O II 



- O VIII obtained from Fig. 03 and assuming that there is 
no interstellar component to these lines. The spectra from 
9th June and 24th November both have sufficient signal 
to noise ratios to provide useful limits. These limits are 
given in Tabled along with the limits on column density 
assuming a velocity width a > 100 kms" 1 . The inferred 
total column density for O II - O VIII is < 5.2 x 10 17 cm" 2 
on the 9th June and < 2.5 x 10 17 cm~ 2 on the 24th Novem- 
ber. 



3.2 Continuum shape 

Having established that the soft X-ray spectrum of 
LMC X-3 is essentially free of discrete emission and ab- 
sorption features except for the interstellar absorption, we 
now consider the shape of the continuum. For all the fits 
performed in this section we include a cold absorbing col- 
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Figure 3. The RGS spectra of LMCX3 from the five differ- 
ent observations. To improve presentation the four brightest 
spectra have been rebinned by a factor of 8. 



Table 2. 3cr upper limits to equivalent widths of O II — O VIII 
Is — 2p absorption lines. The upper limits to column density 
have been calculated assuming velocity width <r > 100 kms -1 . 
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24 Nov — 


Ion 
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EW 
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EW 
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(mA) 


(10 16 cm~ 2 ) 


(mA) 


(10 16 cm~ 


O VIII 


18.97 


< 8 


< 0.8 


< 17 


< 2.2 


O VII 


21.60 


< 37 


< 11 


< 36 


< 9.0 


O VI 


22.01 


< 34 


< 7.8 


< 20 


< 1.6 


O V 


22.33 


< 36 


< 9.4 


< 30 


< 4.1 


O IV 


22.78 


< 13 


< 2.4 


< 23 


< 5.3 


O III 


23.11 


< 47 


< 21 


< 14 


< 2.9 


O II 


23.30 


< 34 


< 6.3 


< 40 


< 8.7 
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Figure 9. Equivalent widths of emission and absorption lines 
of up to 1000 km s — 1 FWHM from the line search described in 
Section 13. 1.21 negative equivalenth width refers to absorption 
lines. The black curve shows the best fit equivalent width as a 
function of wavelength; the grey curves show the 3ar (Ax 2 = 9) 
upper and lower limits. 
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umn fixed at our measured value of 3.8 x 10 20 cm -2 , and 
all the results are tabulated in Table [3] 

We started by fitting a disk-blackbody model to our 
spectra. Only the very faint spectrum from 19th April is 
acceptably fitted by this model, which falls off too rapidly 
below 10A to fit the other spectra. At shorter wavelengths, 
BHC spectra typically require a powerlaw component in 
addition to (or instead of) a disk blackbody (e.g. Ebi- 
sawa et al. 1996). We therefore added a powerlaw to the 
model, with the photon index constrained to lie within the 
range of photon indices seen at higher energy in LMC X- 
3: 1.5 < T < 4.5 (Wilms et al. 2001, Wu et al. 2001). 
This resulted in statistically acceptable fits to the Febru- 
ary, March and April observations, but not to the higher 
signal to noise June and November spectra. In all cases the 
best fit photon index converged to the minimum of the ac- 
ceptable range, indicating that the powerlaw produces too 
many photons at long wavelength, where the data are al- 
ready quite well reproduced by the disk-blackbody compo- 
nent. Therefore we substituted the powerlaw component 
for Compton scattered emission from a hot corona, mak- 
ing use of the THCOMPDS 1 model (Zycki, Done & Smith 
2001). This model upscatters photons to produce a power- 
law shape; the seed photons for the Compton scattered 
component were assumed to come from the disk black- 
body, and the electron temperature of the scatterer was 
fixed to 100 keV (the RGS energy range is too limited for 
the fit to be sensitive to this parameter). This model pro- 
duced better fits to all the spectra, and all of the fits are 
acceptable at the 3<r confidence level; Fig. 1101 shows RGS 
spectra compared to the best fit model components. 



4 DISCUSSION 

In the preliminary analysis of XMM-Newton EPIC and 
RGS spectra, Wu et al. (2001) argued that mass transfer 
mainly occurs in LMC X-3 by Roche lobe overflow rather 
than via a strong stellar wind because the line of sight 
absorption is small (Nh < 10 21 cm -2 ). With the more 
detailed analysis presented here we can rule out accretion 
from a stellar wind with much higher confidence. 

For a continuous spherical wind leaving a body of 
mass M and radius R at the escape velocity, the hydrogen 
column density Nh is related to the mass-loss rate by: 

^ = 24.7m p VGMR x N H (1) 

where m p is the mass of a proton. This expression can 
be used to obtain a conservative upper limit to the wind 
mass-loss rate, because if we assumed a /3-law velocity pro- 
file (dr/dt oc [1- R„/r] ,r > 1.04R,), as used to model 
the line-driven winds from hot stars (Pauldrach, Puis & 
Kudritzki 1986), instead of a constant velocity wind, the 
mass-loss rate for a given column density would be lower 
than that given in equation Q . In Section 13.1.11 we ob- 
tained an upper limit of 1.4 x 10 20 cm -2 for the hydrogen 
column density intrinsic to LMC X-3, based on the neu- 
tral O edge and assuming Solar metallicity. If we instead 
assume a metal abundance of 0.4 x Solar as expected in 

1 THCOMPDS was provided by Piotr Zycki, and is available 
from http://www.camk.edu.pl/~ptz/relrepr.html 



the LMC (Caputo et al. 1999), we find an upper limit of 
Nh < 3.5 x 10 20 cm -2 . We can also estimate the hydrogen 
column density associated with O II - O VIII using the 
limits obtained in Section [3.1.21 from the 24th November 
spectrum. Assuming 0.4 x Solar abundance, we obtain 
Nh < 7.4 x 10 20 cm -2 . Summing these two upper limits 
to obtain the total hydrogen column density associated 
with O I - O VIII, we find N H < 1.1 x 10 21 cm" 2 . As- 
suming that the companion star has mass M < 7Mq and 
radius R < 5Rq (Soria et al. 2001) and substituting into 
equation Q, we find that dM/dt < 8 x 10 17 g s" 1 . This 
is insufficient to power the X-ray emission except when 
LMC X-3 is in its very weakest state (e.g. April 2000). 

However, we have so far neglected wind material in 
which O is completely ionized. When the luminosity of the 
system is as high as 10 38 erg s , the ionization parame- 
ter (£ = L/nr 2 ) at the base of a homogeneous, constant- 
velocity wind from the secondary star is ~ 3 x 10 25 /Nh 
where L is luminosity in erg s" 1 , n is the ion density 
in cm -3 , and r is the binary separation, ~ 10 12 cm. 
Thus for mass-loss rates < 3 x 10 20 g s _1 we expect 
£ > 100 erg cm s _1 , implying that the majority of O 
in the wind is completely ionized (Kallman & McCray 
1982), and therefore undetectable by X-ray absorption line 
spectroscopy. For a wind with a /3-law velocity profile the 
wind is much denser near the base, but at a distance of 
i?»/4 from the stellar photosphere, when the wind has only 
achieved about a fifth of its terminal velocity (assuming 
P ~ l)i £ reaches ~ 10 25 /Nh- A stellar wind cannot be 
driven by line-pressure at such a high ionization state, be- 
cause the ions that provide most of the acceleration (Ab- 
bot 1982) will not be present (McCray & Hatchett 1975). 
On the other hand, we can rule out a mass-loss rate of 
> 3 x 10 20 g s _1 , because the column density would be 
high enough, (and £ low enough) that absorption lines 
from O would be very apparent in the RGS spectra. Thus 
the lack of absorption lines in the RGS spectra implies 
that a line-driven stellar wind is not the primary source of 
accretion material for LMC X-3 except, possibly, when it 
is in a faint state. By comparing the variability properties 
of LMC X-3 to LMC X-l, Nowak et al. (2001) also con- 
clude that LMC X-3 is not accreting from a stellar wind. 
Roche lobe overflow is therefore the most likely mode of 
mass transfer in LMC X-3. 

In order to explain the state changes observed in 
LMC X-3 along with the long term ~ 100 day periodic 
variability, Wilms et al. (2001) propose that the accretion 
rate could be controlled by an accretion disk wind-driven 
limit cycle as discussed by Shields et al. (1986). The op- 
tical photometric variability presented by Brocksopp et 
al. (2001) is also consistent with this hypothesis. In this 
scenario, the central X-ray source in LMC X-3 drives a 
Compton heated wind from the accretion disk that is suf- 
ficiently strong that it interrupts its flow of fuel. For such 
an instability to take place, the mass loss in the wind 
must exceed the accretion rate by a considerable factor 
(e.g. > 30, Shields et al. 1986) . At distances from the cen- 
tral source greater than the outer radius of the accretion 
disk, the Compton heated wind is approximately spherical 
(Begelman & McKee 1983) and so we can use equation 
to relate the mass loss rate of an accretion disk wind 
to the column density. Taking the radius of the accretion 
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Figure 10. Best fit disk-blackbody + Comptonised disk models to the RGS spectra. The dashed black line shows the disk-blackbody 
component, and the dotted black line shows the Comptonised component; the solid line is the sum of the two. The data are shown 
in grey. The data/model ratio is shown in the lower panels; for clarity the error bars are not shown on the data/model ratio plots 
except for the 19 April spectrum. 



disk to be ~ 10 12 cm, and the mass of the primary to be 
~ WMq (Soria et al. 2001), we find that an accretion disk 
wind capable of driving the long-term variability will have 
Nh > 3 x 10 22 cm -2 when the luminosity of the source 
is ~ 10 38 erg s _1 . If we assume a wind temperature of 



~ 10 s K (Begelman, McKee & Shields 1983), more than 
99.9% of O is in the form of O IX, and the column density 
of O VIII can therefore be < 10 16 cm -2 , consistent with 
the upper limits given in Table H Thus we cannot rule out 
that the periodic variability in LMC X-3 is driven by an 
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Table 3. Model fit results. Errors are quoted at 95% confidence for 1 interesting parameter. In cases for which a parameter has 
reached a limit on the allowed range of values while A% 2 < 4 the limit has been marked with an '*'. Fluxes are given in the 0.3 - 2 
keV RGS energy band 

Disk blackbody 



Observation 



kT in 
(keV) 



(10 



Flux (disk) 
-io er g g — i cm - 



Prob 



07 Feb 

08 Mar 
19 Apr 

09 Jun 
24 Nov 



0.63 
0.58 
0.41 
0.34 
1.04 



+0.02 
0.02 
+0.01 
0.01 
+0.15 
0.07 
+0.01 
0.01 
+0.01 
0.01 



2.09 
1.99 
0.007 
0.68 
6.72 



+0.16 
0.13 
+0.09 
0.10 
+0.010 
0.005 
+0.02 
0.02 
+0.24 
0.35 



1169/986 
1161/997 

15/16 
2872/998 
1344/1001 



4.5 x 10~ 5 
2.4 x 10" 4 
0.53 


1.9 x 10~ 12 



Disk blackbody + powerlaw 



Observation 



kT in 
(keV) 



(10 



Flux (disk) 
-io er g g — i cr 



Flux (powerlaw) 



(IO" 



erg s cm 



Prob 



07 Feb 

08 Mar 
19 Apr 

09 Jun 
24 Nov 



67+ 05 
UiD '-0.02 

59+ 01 

0- 28±s;ii: 

0.281^1 

1- i8±S:Si 



1.25" 1 



-+0.30 
-0.27 

1.30±°;» 
0.005™ 



0.28 
6.04 



0.005* 
+0.01 
0.01 
+0.31 
0.29 



-,+0.12 

-0.00* 
-,+0.05 

-0.00* 
-,+3.50* 
-0.00* 
-,+0.02 
-0.00* 

1 K0+ - 09 
1 - ou -0.00* 



1.50: 

1.50j 
1.50j 
1.50"* 



0.89 
0.73 
0.002 
0.43 
0.73 



+0.14 

0.20 
+0.11 
0.08 
+0.004 
-0.002* 
+0.02 
0.03 
+0.08 
0.14 



1045/984 
954/995 

13/14 
1132/996 
1210/999 



0.085 
0.82 
0.54 
1.7 x 10" 
4.7 x 10- 



Disk blackbody + Comptonised disk (thcompds) 



Observation 



kT in 
(keV) 



(10 



Flux (disk) 
" 10 erg s — 1 cm- 



Flux (Comptonised) X 2 l v 

(10~ 10 erg s" 1 cm -2 ) 



Prob 



07 Feb 

08 Mar 
19 Apr 

09 Jun 
24 Nov 



n o 9 +0.05 
u, ° -0.07 
+0.03 
0.03 
+0.24 
-0.14* 

0.26+ - 01 



0.23 



+0.14 



0.36 
0.24 



0.55 



0.01 
+0.02 
0.01 



,+0.07 
-0.42 
003+ 006 

0.43 ± 0.02 

+0.1 
0.0* 



o.6r 



o.o** 



1.60 
1.50 
1.50 
1.50 
1.50 



+0.38 

0.10* 
+0.42 

0.00* 
+3.50* 
0.00* 
+0.33 

0.00* 
+0.03 
0.00* 



1.92 
1.42 
0.005 
0.28 
6.76 



+0.29 
0.29 
+0.49 
0.13 
+0.003 
-0.005* 
+0.04 
0.01 
+0.16 
0.16 



1024/984 

902/995 

12/14 
1079/996 
1127/999 



0.18 
0.98 
0.59 
0.034 
2.8 x 10- 



accretion disk wind. However the column density in any 
accretion disk wind must be small enough that that it is 
transparent to scattering: thermal emission from a wind 
with column density > 10 24 cm -2 would dominate the X- 
ray luminosity of the system, giving rise to prominent Ka 
lines from Fe XXV and Fe XXVI. This is not observed 
in LMC X-3 (Wu et al. 2001). Thus even if there is an 
accretion disk wind, it imparts no significant features or 
distortions on the X-ray spectrum of LMC X-3. 

The resolution and throughput of the RGS coupled 
with the lack of X-ray absorption in LMC X-3 provide 
a clearer picture of the soft X-ray continuum spectral 
shape than has ever been possible before. Excluding the 
very faint April spectrum, the spectral fitting performed 
in Section [3.21 demonstrates that a single disk blackbody 
is unable to describe the 0.3-2 keV spectral shape even 
when LMC X-3 is in the high/soft state, mainly because 
the model falls off too fast at short wavelength. Adding 
a powerlaw component results in much improved (though 
for the 24th November still unacceptable) fits. However, 
the fitted powerlaw slopes consistently bottom out at the 
hardest allowed values in the fitting (r = 1.5) and are con- 
siderably harder than the powerlaw slopes found over the 



0.3-12 keV energy range with XMM-Newton EPIC (Wu 
et al. 2001), even when the EPIC and RGS data were ob- 
tained simultaneously. 

A visual inspection of Fig. |3] offers some insight into 
the reason for the discrepant powerlaw parameters: long- 
ward of 23 A the spectra have almost the same shape, 
running approximately straight and parallel as expected 
from a disk blackbody with a varying inner temperature. 
However, a powerlaw, even with a photon index as hard as 
r = 2, has a softer spectral shape than a disk blackbody 
longward of the peak. Thus for a powerlaw with T > 2 to 
make a significant contribution at A < 10A it would have 
to make a larger contribution at longer wavelength, and 
therefore the spectra should deviate significantly from the 
disk blackbody shape at A > 23A. Therefore the straight 
parallel spectra at A > 23A, and the low best-fit photon 
indices, imply that the component which is observed to 
have a powerlaw shape at < 6 A (i.e. at E > 2 keV) does 
not continue with the same powerlaw shape through the 
soft X-ray band. 

The powerlaw component is normally assumed to be 
the Compton upscattered emission from a hot corona 
(Sunyaev & Titarchuk 1980). The Compton scattered 
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spectrum is expected to deviate from a powerlaw shape 
close to the energy of the seed photons. Thus the absence 
of the powerlaw component in the soft X-ray is consis- 
tent with the seed photons for the Comptonisation com- 
ing from the disk blackbody component. This is predicted 
by successful models for black hole binary spectral states 
(e.g. Haardt & Maraschi 1991, Haardt & Maraschi 1993, 
Chakrabarti & Titarchuk 1995, Haardt et al. 1997), in 
which the temperature and optical depth of the corona de- 
pend on the amount of cooling via Compton scattering off 
soft X-ray photons from the disk. This result is borne out 
by the improved fits which are obtained when the power- 
law is exchanged for a Comptonised disk blackbody. How- 
ever, even with this model, the photon index which de- 
scribes the Comptonised spectrum towards shorter wave- 
lengths is lower than the photon indices obtained over 
a larger energy range using EPIC. Furthermore, for the 
brightest spectra the soft X-ray flux from the Comptonised 
component dominates the flux directly emerging from the 
disk. This suggests some further complexity in the spectral 
shape at the short wavelength end of the spectra (the disk 
blackbody and Comptonised components have a similar 
spectral shape at the long wavelength end). One possi- 
ble explanation could be that in the soft state, the soft 
component in LMC X-3 falls off slower than a disk black- 
body towards short wavelengths due to electon scattering 
near the surface of the disk. Such a broadening of the soft 
component spectral shape has been observed in the soft 
X-ray transients (e.g. GS 2000+25, GRO 1655-40) stud- 
ied by Zycki et al. (2001) using Ginga and RXTE spectra. 
However, the RGS does not extend to short enough wave- 
lengths for a proper investigation of this, which we defer 
to a future study using the EPIC instruments which cover 
a wider energy range. 



5 CONCLUSIONS 

We have presented XMM-Newton RGS spectra of LMC X- 
3 covering the full range of spectral states observed in this 
source (more than 3 orders of magnitude in RGS count 
rate). The spectra are completely dominated by contin- 
uum emission, and are devoid of notable spectral features 
except for the interstellar O I edge. The depth of the edge 
is consistent between the different observations, and is 
consistent with the Galactic interstellar column density 
of 3.2 x 10 20 cm~ 2 in the direction of LMC X-3. The upper 
limit to the edge depth implies Nh < 4 x 10 20 cm -2 asso- 
ciated with neutral O, intrinsic to LMC X-3. Upper limits 
on the equivalent widths of Is — 2p resonant absorption 
lines from O II - O VIII imply N H < 7.4 x 10 20 cm- 2 for 
gas associated with partially ionized O. This implies that 
a line-driven stellar wind off the companion star can not 
supply a significant fraction of the accreting material ex- 
cept during the very low luminosity states. Therefore the 
majority of accretion in LMC X-3 probably takes place 
through Roche Lobe overflow. 

A multi temperature disk blackbody is a good repre- 
sentation of the soft X-ray spectrum of LMC X-3 longward 
of the disk blackbody peak. Thus the powerlaw compo- 
nent which dominates at higher energies does not continue 
through the soft X-ray range with a powerlaw shape. This 



suggests that the disk blackbody component supplies the 
seed photons which are Compton upscattered to form the 
powerlaw at higher energies, as in conventional models for 
BHC spectral states. 
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